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ABSTRACT* 

The  interrelation  between  field  programs  and  atmospheric  modeling 
is  discussed  in  the  context  of  studies  of  weather  modification.  In 
particular,  numerical  models  of  convective  clouds  are  considered.  It 
is  shov  liow  activities  in  both  field  programs  and  modeling  affect 
decisions  concerning  each  other,  culminating  in  the  use  of  models  for 
day-to-day  "go"  and  "no  go"  decisions  concerning  the  field  programs. 
Illustrations  are  given  from  work  done  at  RAND  and  elsewhere.  A  few 
typical  results  from  the  RAND  cumulus  dynamics  model  are  presented. 

This  talk  was  prepared  for  the  Sixth  Skywater  Conference,  spon¬ 
sored  by  the  Bureau  of  Reclamation  in  Denver,  10-11  February  1970. 
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should  not  be  interpreted  as  reflecting  the  views  of  The  RAND  Corpora¬ 
tion  or  the  official  opinion  or  policy  of  a n;  of  its  governmental  or 
private  research  sponsors.  Papers  are  reproduced  by  The  RAND  Corpora¬ 
tion  as  a  courtssv  to  msibars  of  its  staff. 
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ATMQSPHERIC  MODELING,  FIELD  PROGRAMS,  AND  DECISION  SYSTEMS 

I 

It  was  something  of  a  surprise  to  me  to  find  myself  scheduled  to 
talk  on  the  role  of  modeling  and  decision  systems,  for  that  appears  to 
involve  close  coupling  of  field  programs  with  modeling,  and  1  have  had 
no  direct  connection  with  a  field  program.  But,  as  I  hope  to  make 
clear,  there  is  always  such  a  connection,  whether  direct  or  indirect, 
and  it  is  of  utmost  importance  to  both  activitiea.  The  modeler  must 
be  strongly  concerned  with  field  programs  even  if  he  does  not  engage 
in  them  himself.  I  am  glad  that  our  moderator  said  that  there  was 
little  distinction  between  the  topics  for  yesterday  afternoon  and  this 
morning,  for  my  subject  is  the  interrelation  between  field  programs 
and  modeling  as  it  ia  expressed  in  decision  making. 

We  have  already  heard  many  good  reasons  for  engaging  in  modeling 
along  with  field  programs.  One  of  those  reasons  is  to  use  models  to 
help  make  decisions  concerning  the  field  programs.  In  its  most  evi¬ 
dent  application  thla  can  mean  using  a  model  the  way  the  Experimental 

1  2 

Meteorology  Laboratory  does  in  Miami,  '  or  the  way  the  MRI  group  does 
3 

in  Arizona  —  to  determine  which  clouds  are  the  best  for  seeding  or 
which  days  are  the  most  favorable.  But  this  usage  may  come  rather 
late  in  the  scheme  of  things.  Actually,  decisions  must  be  made  all 
along  the  line,  both  on  the  modeling  and  on  the  field-program  sides, 
and  simultaneous  activity  on  both  aides  offers  the  best  hope  of  success. 

Both  field  programs  and  modeling  may  be  divided  into  "development" 
and  "production"  stages,  although  those  words  may  have  somewhat  dif¬ 
ferent  connotations  in  different  contexts.  In  the  context  of  the 
present  meeting,  the  "development"  stage  of  a  field  program  could  in¬ 
clude  observations  of  clouds  to  learn  more  about  the  processes  going 
on  In  them,  end  it  could  also  inelude  experimental  attempts  at  modi¬ 
fication.  Later,  when  we  have  developed  a  reasonably  good  knowledge 
of  the  results  of  certain  modification  activities,  we  can  move  to  the 
"production"  stage  and  engage  in  these  activities  for  practical  rea¬ 
sons,  such  as  increasing  rainfall  or  snowfall  in  a  particular  region. 

In  some  programs,  the  "production"  stage  could  come  earlier. 
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in  the  case  of  modeling  (and  I  refer  only  to  numerical  modeling, 
though  many  of  the  same  things  could  be  said  for  laboratory  modeling), 
the  "development"  stage  covers  the  initial  formulation  and  program¬ 
ming.  It  must  also  include  some  type  of  check  to  determine  what  re¬ 
lation  the  results  bear  to  reality.  Here  ;  the  first  obvious  spot 
where  a  field  program  can  exercise  a  decision-making  influence  on  a 
model.,  The  checking  of  model  results  against  actual  observations 
tells  us  whether  we  are  on  the  right  track  in  our  model  building  and 
suggests  changes  to  be  made  in  the  model.  I  will  have  more  to  say 
about  this  kind  of  checking  later.  At  thia  stage,  declalon-maklng 
can  go  in  the  opposite  direction,  too,  for  the  characteristics  of  the 
model  under  consideration  may  suggest  the  types  of  observations  that 
are  most  useful.  Frequently  in  developing  a  model  one  is  confronted 
with  fne  need  to  know  the  value  of  e  quantity  that  is  not  or  cannot 
be  observed.  Sometimes  theory  will  suggest  an  obaervabla  substitute 
quantity,  but  at  other  times  the  field  experimenter  is  spurred  to 
develop  new  equipment  or  new  uses  for  old  equipmant  to  satisfy  the 
modeler. 

Only  after  thorough  checking  of  the  model  results  egainet  actual 
observations  can  the  "production"  stage  of  a  model  be  entered.  One 
aspect  of  this  stage  la  the  operational  use  of  a  modal  for  day-to-day 
implementation  of  the  field  program,  as  mentioned  previously.  Another 
aspect  is  the  use  of  the  model  itself  in  the  direct  study  of  the  thing 
being  modeled. 

Atmospheric  processes  ere  extremely  complex,  and  ona  of  tha  func¬ 
tions  of  both  field  experimentation  and  modallng  la  to  rendar  tha  com 
plexlty  more  intelligible.  Many  of  the  proceeees  not  well  under¬ 
stood,  which  makes  the  design  of  the  model  or  experiment  difficult, 
but  frequently  a  we  11  -designed  model  wij. ’  lead  us  to  an  understanding 
of  such  processes.  The  particular  advantage  of  e  model  in  this  re¬ 
spect  is  its  ability  to  isolate  tha  effects  of  individual  proceaaaa 
or  parameters.  This  was  discussad  yeatarday  by  Orville.  It  is  easy 
to  run  through  a  serins  of  numerical  experiments,  varying  one  compo¬ 
nent  of  the  model  to  see  how  this  affects  the  total  results.  However, 
if  this  ie  to  be  more  than  a  mathamatical  exercise,  it  must  be  tied 


to  observations  of  the  real  atmosphere.  Hence,  modeling  must  not  pro¬ 
ceed  independent  of  field  experimentation.  On  the  other  hand,  a  field 
program  not  based  on  the  understanding  of  theory  that  comes  from  model¬ 
ing  could  flounder  and  produce  great  masses  of  undigestible  data.  The 
implicitness  of  atmospheric  processes  makes  it  virtually  impossible 
for  a  field  program  to  study  one  process  in  isolation;  hence  the  re¬ 
course  to  models. 

Thus,  we  see  that  models  and  field  programs  should  each  contribute 
to  decisions  concerning  the  other  in  their  earliest  stages  of  develop¬ 
ment.  Throughout  the  process,  decisions  concerning  the  direction  of 
work  on  a  model  should  be  based  on  field  observations  as  much  as  on 
theory.  Once  a  model  achieves  a  fair  degree  of  realism,  even  if  only 
in  one  respect,  such  as  the  ultimate  height  of  convection,  it  can  be 
used  as  a  basis  for  routine  day-to-day  decisions  about  Held  opera¬ 
tions.  At  a  higher  degree  of  development,  the  model  can  be  used  to 
study  the  atmosphere  dlreetly,  but  still  with  continual  verification 
from  field  observation. 
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II 


In  Che  many  interactions  that  I  have  touched  upon,  the  flow  of 
information  betwecr.  the  experimenter  and  the  theoretician  may  be  ex¬ 
tremely  complex.  When  both  are  represented  in  a  combined  operation, 
matters  are  simplified.  A  good  example  of  such  a  combined  operation 
is  the  Experimental  Meteorology  Laboratory.  The  flow  of  information 
may,  however,  be  largely  through  the  literature.  A  good  example  of 
this  is  the  parameterlsation  of  cloud  microphysics  by  Kessler.  There 
is  a  great  deal  of  knowledge  available  concerning  the  theory  of  micro- 
physical  processes  in  clouds,  and  several  numerical  models  have  been 
developed  to  describe  the  changes  in  the  drop-sise  spectrum.  Concur¬ 
rently,  models  have  been  constructed  to  describe  the  dynamics  of  whole 
clouds,  but  because  a  gross  mismatch  of  both  time  and  spare  scale, 
these  types  of  models  have  not  yet  successfully  been  combined.  What 
was  needed  was  a  parameterlsation  whereby  the  microscale  processes 
could  be  described  with  reasonable  accuracy  in  terms  of  cloud-scale 
parameters.  Making  full  use  of  observations  from  field  and  laboratory 
experiments  quoted  in  the  literature,  together  with  physical  theory, 

Kaesler  produced  the  required  parameterlsation  for  three  of  the  most 
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important  mlcrophyolcal  processes.  '  Upon  its  publication,  this 
parameterlsation  was  picked  up  by  many  cloud  modelers  and  Incorporated 
into  their  own  models,  often  with  variations,  such  as  the  alternate 
expression  for  conversion  proposed  by  Berry. ®  The  incorporation  of 
this  parameterlsation  has  so  Improved  the  models  that  most  of  them  now 
contain  it  or  something  very  much  like  it.  That  the  models  were  im¬ 
proved  has  been  verified  by  field  observation. 


-6- 


III 

The  ultimate  use  to  which  a  model  will  be  put  determines  to  a 
large  extent  what  kind  of  model  it  will  be.  If  the  model  is  to  be  a 
general  reoearch  tool,  designed  to  reveal  new  information  about  the 
workings  of  the  atmosphere,  it  must  be  as  comprehensive  as  possible. 
Ideally  this  would  mean  going  back  to  the  basic  equations  of  hydro¬ 
dynamics,  thermodynamics,  physics,  etc.,  to  describe  all  possible 
processes,  and  solving  these  equations  without  approximation.  This 
is  an  unattainable  goal  on  several  counts.  Compromises  must  be  made 
all  along  the  line.  However,  it  la  possible  to  specialize  the  equa¬ 
tions  to  a  particular  scale  of  motion  (thereby  simplifying  them  sig¬ 
nificantly),  parametcrlse  the  interaction  with  other  scales,  and  solve 
them  by  numerical  methods.  Even  with  drastic  pruning,  such  a  model  is 
usually  found  to  tax  the  capabilities  of  the  largest  and  most  advanced 
computers  permitted  by  the  existing  technology.  Cloud  models  fitting 
this  description  have  been  <n  use  now  for  nearly  ten  years,  and  their 
requirements  have  always  kept  pace  with,  or  remained  slightly  ahead 
of,  whatever  computational  facilities  have  been  available.  As  the 
computers  grew,  so  did  the  models.  Perhaps  now  is  the  time  to  adopt 
some  radically  new  approach,  but  I  cannot  discuss  this  today. 

Throughout  all  tnis  growth  it  has  been  necessary  to  .is  the  models 
to  reality  by  extensive  comparison  of  their  results  with  field  obser¬ 
vations.  The  model  discussed  by  Orville  is  of  this  type,  and  he  has 
the  field  programs  of  the  South  Dakota  School  of  Mines  to  work  with.7 
The  model  I  developed  is  also  of  this  type,  and  I  have  had  sooe  access 

to  the  results  of  the  field  programs  in  the  Caribbean  of  the  Naval  Re- 
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search  Laboratory  and  the  Experimental  Meteorology  Laboratory.  *  Al¬ 
though  these  sad  other  similar  models  have  aided  in  the  understanding 
of  meny  of  the  facets  of  cloud  dynamics,  they  are  not  sufficiently 
comprehensive  to  be  able  to  supplant  field  programs,  nor  are  they  likely 
to  becoaw  so.  Their  use,  however,  can  suggest  certain  types  of  field 
programs  that  might  be  worthwhile  to  undertake  and  others  that  would 
probably  be  useles*. 
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For  some  purposes  a  much  more  simplified  model  can  be  useful. 
Typically,  this  sort  of  model  considers  only  one  dimension  rather  than 
two  or  three,  and  it  specifies  the  cloud  to  be  of  some  idealized  form, 
such  as  a  bubble  or  a  starting  plu'se.  A  large  degree  of  parameteriza¬ 
tion  is  required,  not  only  for  the  microphysics,  but  also  for  processes 
at  the  larger  scale.  Such  models  tend  to  be  deficient  in  dynamics  and 
can  offer  little  toward  understanding  the  circulations  in  and  around 
clouds.  Nevertheless,  they  have  the  virtue  of  fast  running  on  a  small 
computer,  and  they  are  able  to  predict  with  considerable  accuracy  cer¬ 
tain  features  of  cloud  development,  notably  •'he  maximum  height.  These 
models  f, re  well  suited  to  direct  use  in  conjunction  with  a  field  pro¬ 
gram,  as  has  been  well  demonstrated  by  the  Experimental  Meteorology 
Laboratory  and  by  MR1. 
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IV 


It  seems  to  be  de  vigour  to  talk  a  little  about  one's  own  pro¬ 
gram,  so  I  will  do  so.  I  will  not  attempt  to  describe  it  in  detail 
or  give  a  lot  of  results,  but  will  indicate  how  we  try  to  keep  it 
tied  to  reality. 

At  RAND  we  have  no  direct  connection  with  a  field  program,  and 
we  have  elected  to  concentrate  on  the  large  hydrodynamical  models. 

But,  as  I  have  mentioned,  such  models  must  continually  be  compared 
with  field  observations  to  assure  that  they  are  doing  what  they  are 
intended  to  do.  One  of  the  means  ve  have  adopted  is  to  use  only  real 
atmospheric  soundings  for  initial  data.  This  somewhat  complicates  the 
results,  for  an  idealized  initial  condition,  such  as  a  constant  lapse 
rate,  for  example,  would  lead  to  simpler  patterns  for  all  the  vari¬ 
ables.  But  we  think  that  the  closer  correspondence  to  what  is  ob¬ 
served  in  nature  is  worth  the  added  complication.  In  choosing  the 
initial  sounding  we  attempt  to  find  conditions  for  which  verifying 
data  are  available.  This  has  generally  meant  using  Caribbean  sound¬ 
ings  for  days  on  which  tue  instrumented  aircraft  of  the  Naval  Research 
Laboratory  or  the  Experimental  Meteorology  Laboratory  were  active. 

Even  so,  the  verifying  data  can  never  be  as  extensive  as  we  would  like. 
The  best  we  can  usually  hope  for  is  a  few  penetrations  showing  liquid 
water  content  and  temperature,  together  with  height  of  base  and  top 
of  the  cloud  and  cloud  width.  The  age  of  the  cloud  being  observed  is 
seldom  known,  so  the  choice  of  a  time  step  in  the  simulation  Is  some¬ 
what  arbitrary.  Temperatures  and  dew  points  are  hard  to  compare,  for 
the  model  works  with  departures  from  the  basic  state,  and  that  basic 
state  cannot  wall  be  observed  in  practice.  A  most  unfortunate  circum¬ 
stance  le  that  the  aircraft  are  not  well  able  to  measure  horizontal 
and  vertical  components  of  wind  throughout  the  cloud. 

Nevertheless,  we  have  been  able  to  get  a  fair  amount  of  data  for 
comparison.  As  an  illustration,  Figure  1  shows  a  comparison  of  compu¬ 
tation  with  four  clouds  observed  by  NRL  aircraft  near  Puerto  Rico  on 
one  <i»y.  As  1  have  mentioned,  the  correspondence  of  temperature  de¬ 
part  ire  and  dew  point  is  not  good,  but  this  is  related  to  the  diffi- 
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loud 


Property 

A 

— 

B 

C 

D 

Notes 

Height  of  penetration  (ft! 

6  ,600 

7 ,100 

7,350 

7,700 

Temperature  departure  (°C) 

-0.A 

-1.0 

-1.2 

-2.2 

a 

(inside  temperature  minus 

1.8 

-1  . 1 

0.0 

-1.0 

fc 

outside  temperature) 

1.5 

1.5 

1.3 

1.2 

c 

Evaporator  dew  point  (°C) 

13.3 

13.5 

12.8 

10.6 

a 

16.9 

15.8 

13.8 

9.0 

b 

16.5 

16.0 

15.8 

15.5 

c 

-3 

Liquid  water  cont  nt  (g  m  ) 

2.0 

2.3 

2.3 

1.2 

d 

1.4 

2.2 

0.0 

1.2 

c 

2.1 

2.1 

1.5 

0.6 

b 

2.0 

2.4 

2.6 

3.0 

c 

Width  of  cloud  (ra) 

1,200 

600 

1  ,000 

900 

a 

1 ,070 

1,040 

430 

120 

b 

610 

790 

910 

910 

c 

Height  of  c  leud  base  (ft) 

2  ,800 

2  ,800 

2  ,800 

2,800 

a 

3,200 

3  ,200 

3,200 

3,200 

b 

3  ,200 

3,200 

3,200 

3,200 

c 

Height  of  cloud  top  (ft) 

-10,000 

-10,000 

-10,000 

-10,000 

a 

7  ,900 

7,900 

7  ,900 

7  ,900 

b 

9,100 

9,100 

9,100 

9.100 

c 

a  Aircraft  observation, 
b  Computation  at  20  minutes, 
c  Computation  at  23  minutes, 
d  Observed  by  NRL  total  water  instrument, 
e  Observed  by  Johnton-Wi 1 1 i ams  liquid-water  meter. 

Fig.  1  --  Cbtervad  and  computed  value*  of  cloud  properties. 
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culties  In  relating  the  ages  of  the  real  and  simulated  clouds,  lack 
of  a  proper  base  line,  and  lack  of  knowledge  of  the  updrafts  and 
downdrafts  in  the  real  clouds.  I  conclude  from  this  that  improve¬ 
ments  are  desirable  both  in  the  numerical  model  and  in  the  observing 
equipment  and  technl'p’es .  The  liquid  water  content  shows  reasonably 
good  agreement  between  computation  and  measurement.  Except  f~r 
Cloud  D  the  measurements  made  with  the  evaporation  dew  point  equip¬ 
ment  agree  with  the  computations  better  than  those  made  with  the 
Johnson-Williams  instrument.  This  is  partly  due  to  their  different 
responses  depending  o  drop  size.  But  careful  examination  of  the 
figures  will  also  show  a  deficiency  in  the  model  related  to  the  fact 
that  this  version  had  no  fall-out  mechanism  for  liquid  water,  aud  so 
tended  toward  large  values  at  the  top  of  the  cloud  and  small  values 
at  the  bottom.  The  width  of  the  cloud,  which  is  such  an  important 
basic  input  to  the  one-dimens ionai  models,  is  of  less  significance 
here.  It  is  somewhat  dependent  on  the  initial  impulse  specified, 
and  its  effect  on  the  ultimate  growth  of  the  cloud  has  not  been  well 
established  in  this  model.  The  height  of  the  cloud  base  in  the  model 
is  determined  both  by  the  initial  sounding  and  the  impulse.  Hence, 
the  experimenter  has  only  partial  control  over  it.  The  ultimate 
height  of  the  cloud  top  can  be  controlled  in  the  same  way,  but  re¬ 
cent  evidence  indicates  that  the  initial  impulse  has  a  rather  strong 
effect. 

A  more  recent  version  of  the  RAND  model  Incorporates  Kessler's 
parameterization  of  microphysics.  In  this  version  liquid  water  is 
divided  into  two  categories:  cloud  water,  which  consists  of  drops 
small  enough  to  move  with  the  air,  and  hydromete  water,  which 
consists  of  drops  large  enough  to  fall  relative  to  the  air.  Figure  2 
shows  how  the  distribution  of  the  two  forms  of  liquid  water  varies 
on  the  central  axis  with  time.  From  the  inception  of  condensation 
out  to  about  11  minutes  only  cloud  water  exists.  At  that  time  the 
liquid  cortent  becomes  great  enough  to  start  the  conversion  to  hydro¬ 
meteor  weter,  and  both  forms  coexist.  Around  15  minutes  rain  starts 
to  fall  from  the  base  of  the  cloud,  first  reaching  the  ground  at 
22  minutes.  Later  the  base  of  the  cloud  lifts,  and  by  the  time  the 


Horizontal  hatching  represents  cloud  water;  vertical  hatching 
represents  hydrometeor  water. 
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cloud  top  reaches  its  maximum  height,  the  small  drops  are  found  only 
In  the  upper  kilometer  or  two.  After  45  minutes,  only  hydrometeor 
water  exists,  and  all  of  it  falls  to  the  ground  before  60  minutes. 

The  availability  of  quantitative  values  of  hydrometeor  water 
suggests  the  computation  of  a  variable  that  has  little  theoretical 
significance  in  itself,  but  is  most  important  in  an  observational 
program;  namely,  radar  reflectivity.  Figure  3  shows  (in  the  dashed 
line)  the  magnitude  of  the  maximum  value  of  radar  reflectivity.  The 
nominal  threshold  for  detectability  is  first  reached  at  15  minutes, 
after  which  the  value  grows  until  34  minutes  and  then  decreases  to 
the  threshold  at  59  minutes.  This  corresponds  with  the  variation  in 
total  liquid  water,  shown  in  the  previous  figure.  The  height  at  which 
the  radar  reflectivity  reached  its  maximum  is  shown  by  the  solid  line, 
and  the  vertical  range  over  which  the  radar  reflectivity  is  large 
enough  for  nominal  detectability  is  indicated  by  shading.  Fortui¬ 
tously,  Saunders  has  published  a  study  based  on  observations  of  r-dar 
reflectivity  over  the  Caribbean  within  a  few  days  of  the  basic  sound¬ 
ing  used  on  the  figure.  The  dotted  line  shows  that  the  value  of  max¬ 
imum  radar  reflectivity  coincided  with  the  computed  value  almost  ex¬ 
actly  for  nearly  ten  minutes,  but  then  leveled  off  more  quickly.  The 
main  difference  is  that  the  observed  height  of  the  maximum  continu¬ 
ously  fell,  whereas  that  of  the  computation  rose  until  just  before  the 
cloud  top  also  reached  its  maximum  and  then  fell.  But  in  the  descend¬ 
ing  stage  the  slopes  of  the  computed  and  observed  curves  were  almost 
identical,  though  they  were  separated  in  time  by  over  ten  minutes. 

The  usual  problem  of  comparison  arises:  there  is  no  wholly  acceptable 
way  of  relating  the  zero  times  of  the  computed  and  observed  clouds, 
but  in  any  case  there  must  br  a  discrepancy  in  this  example.  The 
rates  of  precipitation  growth  match  very  well,  and  so  do  the  rates  of 
descant  of  rain,  but  the  timing  is  off.  We  infer  that  the  radar  echo 
in  the  observed  cloud  first  appeared  when  the  cloud  was  nearing  its 
maximum  vertical  growth,  but  it  appeared  somewhat  earlier  in  the  com¬ 
puted  cloud,  while  there  were  still  substantial  updrafts.  This  may 
be  a  result  of  the  feature  of  Kessler's  parameterization  in  which  the 
rate  of  conversion  from  cloud  to  hydrometeor  water  is  a  linear  function 
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of  cloud  water  content  with  a  threshold.  We  hope  to  try  Berry’s 
formulation  of  the  conversion  process,  and  it  may  lead  to  some  im¬ 
provement  in  this  respect.  Beyond  that  lies  a  long  series  of  im¬ 
provements  in  the  parameterizations  and  other  features  of  the  model, 
always  comparing  the  response  to  any  change  with  observations  of 
real  clouds.  Only  in  that  way  will  model  building  have  any  value 
in  helping  us  to  understand  the  complicated  atmospheric  processes. 
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